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(N ■ 

D ■ ABSTRACT 

^ ■ Here we discuss evolution and broad-band emission of compact (< kpc) lobes in 

young radio sources. We propose a simple dynamical description for these objects, 
• consisting of a relativistic jet propagating into a uniform gaseous medium in the cen- 

tral parts of an elliptical host. In the framework of the proposed model, we follow 
2 i the evolution of ultrarelativistic electrons injected from a terminal hotspot of a jet to 

"to ', expanding lobes, taking into account their adiabatic energy losses as well as radiative 

^ ' cooling. This allows us to discuss the broad-band lobe emission of young radio sources. 

In particular, we argue that the observed spectral turnover in the radio synchrotron 
^ ■ spectra of these objects cannot originate from the synchrotron self- absorption process 

o 

psj , interstellar medium engulfed by the expanding lobes and photoionized by active cen- 



but is most likely due to free-free absorption effects connected with neutral clouds of 



ters. We also find a relatively strong and complex high-energy emission component 
I produced by inverse-Compton up-scattering of various surrounding photon fields by the 

I lobes' electrons. We argue that such high energy radiation is strong enough to account 

^ ' for several observed properties of GHz-peaked-spectrum (GPS) radio galaxies at UV 

and X-ray frequencies. In addition, this emission is expected to extend up to GeV (or 
^ ■ possibly even TeV) photon energies and can thus be probed by several modern 7-ray 
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instruments. In particular, we suggest that GPS radio galaxies should constitute a 
relatively numerous class of extragalactic sources detected by GLAST. 

Subject headings: galaxies: active — galaxies: jets — acceleration of particles — radi- 
ation mechanism: non-thermal 



1. Introduction 



'GHz-peaked-spectrum' (GPS) objects are powerful radio sources whose spectra are inverted 
{Lu oc u^" with a < 0) below peak (or turnover) frequencies fp ~ 0.5 — 10 GHz and whose linear 
sizes are LS < Ikpc. 'Compact steep spectrum' (CSS) objects are similarly powerful inverted- 
spectrum radio sources but with peak frequencies in a lower frequency range when compared to 
the GPS population, Up < 0.5 GHz, and with larger linear sizes, LS ~ 1 — lOkpc. Nuclei of GPS 
and CSS objects can be classified as either radio galaxies, quasars, or Seyfert galaxies of type 1 or 
2. Morphologically, GPS/CSS sources may be reminiscent of a smaller version of classical doubles 
(FR H radio galaxies), with pairs of symmetric lobes present at opposite sides of weak nuclei. In such 
cases, they are called 'compact symmetric objects' (CSOs), if LS < Ikpc, or 'medium symmetric 
objects' (MSOs), if LS ~ 1 — lOkpc. Quite often, however, GPS/CSS sources are characterized 
rather by a 'core-jet' morphology with asymmetric lobes (if present at all). In such cases, it is not 
clear whether they should be classified as 'true' GPS/CSS-es or rather as 'regular' (i.e. extended) 
radio-loud active galactic nuclei (AGNs) viewed in projection. About 10% of radio sources found 
in high-frequency radio surveys belong to the GPS class whereas 30% are classified as CSS objects. 
An extensive review on this issue was presented by lO 'Deal (| 19981 ). 



As shown bv iQ'Dea &: BaumI (|l997l ). a relatively tight correlation between the turnover fre- 
quency and the source's linear size, namely log (t'p/ GHz) = — 0.21 (±0.05)— 0.65(±0. 05) x log (LS'/kpc), 
holds for the investigated parameter range fp = 0.05 — 20 GHz and LS = 0.01 — 20kpc. This uni- 
fies the GPS and CSS populations and suggests that they are both manifestations of the same 
physical phenomenon. Moreover, such a continuous distribution up to the observationally lim- 
ited peak frequenci es ^ 10 GHz sugges t s that there may be an unnoticed population of sources 
with fp > 10 GHz (jO'Dea k BaumI 119971 : IXornikoski et al.ll200d ). These sources are called 'high 



frequ ency peakers' (HFP), and several possible candidates were already selected (jPallacasa et al 



2000). However, the majority of the candidates (especially those that are quasar-h osted) possess 



clear core-jet morphology and thus may be not related to the GPS/CSS phenomenon (jOrienti et al 



20061 ). In this context, monitoring studies may help in performing the pr oper classification since 



little variability (< 10%) is expected for the discussed cl ass of objects (jTornikoski et al.l 12001 
Tinti et ahlbood : Ixorniainen et aklboosl : lOrienti et aklboOTl ). 



The 'true' GPS/CSS sources have to be intrinsically very powerful in radio because Doppler 



and projection effects see m to be, in those cases, rather marginal (iFanti et al 



1994 : ISaikia et al.l Il995l ). with the possible exception of GPS quasars (jStanghellini et al 



1991 



a 



Wilkinson 



et al 



2001 
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Stanghellinill2003l ) . In particular, radio powers of the considered sources at 5 GHz always exceed the 



FRI/FR II division, Li .4nHz > 10 ^^ WHz~^ and reach Li.4ghz > lO^^WMz"! in some cases. As 



shown bv lStanghellini et al.l (|l998l ). 10% — 20% of GPS objects possess, in addit ion, faint extended - 



radio emission, with th e famous 0108+388 being the most obvious example (iBaum et al 



199S 



Marecki et al. 



Stanghellini et al.lll990). Such extended radio halos may reach even Mpc scales ( Schoenmakers et al. 



199C 



20031 ) and are believed to represent fossil structures formed in previous epochs 



of the jet activity. This idea is supported by the presence of a GPS-like radio core in the source 
J1247+6723, which is character ized by a classical 'double-double' (i. e., restarting) large - scale radio 
morphology (jSaikia et al.l 120071 ). On the other hand, as argued bv IStanghellini et al.l (|2005l ). the 
extended emission is most often seen in GPS quasars, which are more likely core-jet like than truly 
compact structures and therefore not necessarily represent fossil lobes in all cases. 

Although there is an emerging agreement that GPS/CSS sources — at least those truly compact 
and not simply shortened by the projection effects — are young versions of extended radio galaxies 
and quasars (see §2.1 below), several key questions regarding these objects remain open. They 
concern, for example, (i) the nature of the absorption mechanism responsible for the observed 
inverted radio spectra at low frequencies, (ii) details of the dynamical evolution and interaction 
with the ambient (galactic) medium, and also (iii) the parameters of the central engine like the 
accretion rate, the nuclear obscuration, etc. Clearly, a detailed analysis of the broad-band emission 
from GPS sources, including recent observations in the X-ray photon energy range, may help to 
answer some of these questions. Here we explore the possibility that young radio galaxies may 
be, in addition, sources of relatively intense 7-ray emission and that detection of such radiation 
(or even the establishment of upper limits to it) by instruments like GLAST, AGILE, H.E.S.S., 
MAGIC, or VERITAS in the GeV-TeV photon energy range can help to constrain the physics of 
this class of objects. In particular, in § 2 we propose a simple and updated dynamical description 
of the evolution of GPS sources, which allows us to discuss in § 3 the expected broad-band emission 
of their lobes, including the GeV photon energy range. Final conclusions are given in § 4 of the 
paper. In a subsequent paper, we carefully select from the literature several GPS radio galaxies 
(where the lobe emission is expected to dominate the total radiative outputs, since the underlying 
relativistic jet and accretion disk emissions are likely to be Doppler- hidden and/or obscured) that 
are detected at X-ray energies and analyze their multiwavelength radiation in the framework of the 
presented model. 



Evolution of GPS Sources 



2.1. Present Understanding 



Since GPS/CSS objects are as powerful as clas sical doubles but muc h smaller, they can be 
either young versions of the extended radio sourc es (iPhilips &: Mute]||l982l ) or examples of radio- 
loud AGNs 'frustrated' by the ambient medium (|van Breugel et al.lll984l ). Efficient confinement 
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of an expanding radio structure by a dense galactic environment was proposed to be associated 
in a natural way with the narrow- line region (NLR). We note that the typical parameters — 
temperature, average number density, and filling factor — of the NLR clouds, as observed in many 
powerful radio galaxies, are T ~ 10^ K, tinlr ~ 10^ — lO'^ cm~^, and (p ~ 10~^, respectively. These 
NLR clouds, distributed around galactic nuclei on kpc scales with total masses up to M^lr ~ 
are embedded within hot, X-ray emitting gaseous halos, whose typical densities (inferred 
from X-ray observations of giant ellipticals) are raisM(^ 1 kpc) ~ 0.1 cm~'^. However, the frustration 
scenario requires total masses of cold a mbient gas in a range 10^*^ — 10^^ Mq within the host galaxies 



( De Young 1993 : Carvalho 1994, 19981). Such significantly dens er environments of GPS /CSS sources 
were indeed claimed previously (jGopal-Krishna &: Wiitalll99lh but are not supp orted by the most 



recent multiwavelength studies (although see recently ICarcia-Burillo et al.ll2007l ). 



In the framework of the youth scenario (IPhilips fc Mutejll982l '). the evolution of GPS/CSS 



sources tow a rd ex tended FRIIs was followed by lCarvalhd (|l985l ) and then by a number of authors. 



/lutell 
198,4 



Fanti et al.l (|l995l ) found that, in order to explain the size distribution of radio galaxies in this 
approach, one has to invoke a decrease of the radio power with increasing linear size for young, 
compact objects. Str ong negative luminosity evolution of GPS/CSS sources (advocated also by 
Readhead et al.lll996l ) would imply that most GPS/CSS objects cannot be precursors of the most 
luminous FRII s but only of the low-power radio galaxies, located close to the FRII/FRI division. 
In this respect. iBegelmanI (j 19961 ) presented a simple evolutionary model that successfully accounted 
for many observed features of the GPS/CSS class. In particular, assuming (i) a power-law density 
profile of the ambient medium pisM{f) with j3 ~ 1.5 — 2 (where r is the distance from the 

galactic nucleus), (ii) constant jet power during the source's lifetime, and (iii) self-similar expansion 
of the overpressured lobe (close to energy equipartition) , IBegelmanI (119961 ) obtained an almost 
constant advance velocity of the hotspots, a decrease of the radio power with size Lr cx r~^'^, and 
a size distribution d N/dloR LS oa r~" ^~^^ with m ~ 0.6 (as implied by observations for LS up to 
hundreds of kpc; see lFanti et al.ll200ll . and references therein). 



O'Dea &: BaumI (jl997l ) noted, however, a more complex size distribution of radio sources, con- 
sisting of a plateau dN/dlogLS oc const for 0.3 kpc < LS < 10 kpc and a power-law tail oc LS^'^ 
for LS > 10 kpc. This could imply an overa bundance of compact radio s ources when compared to 
the number of extended ones. As shown by lReynolds &: BegelmanI (|l997l ). such an overabundance 
can be incorporated into the simple self-similar evolutionary model only if jet intermittency is in- 
troduced f l O^yr- long burst of jet ac t ivity recurring every lO^yrs). Alternatively, as discussed by 
Alexandeii (|200d l and ISnellen et al.l (|2000l ). an additional population of short-liv ed radio sources 
that d ie before reaching > 10 kpc scales may explain the observational results of lO'Dea &: Baum 
(|l997l ). which, however, s till needs to be confirm ed within a sample of sources not affected by the 
projection effects (but see lTinti fc: de Zottill2006l ). The issue of self-similarity in the evolution of ra- 
dio galaxies is still in general debated, and non-self-similar evolutionary models for compact sources, 
enriched by some addition al effects like ene r gy /momentum losses of the jets due to interaction s with 



the surrounding medium (jPe YoungI Il997l : iPerucho &: Marti 



2002 



Kawakatu &: Kind 120061 ). were 
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discussed. Self-similar scenarios were also explored in more detail, enriched by (more consistent 
with observations) King-type ambient medium density profiles pisuir) oc (f -|- (r/rc)^)"^/^ instead 
of a single power-law considered earlier. The presence of a plateau in the ambient gas density within 
the core radius Tc ~ 1 kpc implies that, in the initial state, the radio luminosity of a single s ource 
may even increase with increasing size for LS < and then decrease for LS > rc (|Alexandeiil200Qi : 



Snellen et al.ll2nnm . 



2.2. A Simple Dynamical Model 



All the models describing the evo lution of GPS/CSS sources start from the set of equations 
discussed by iBegelman &: CiofBl (j 19891 ) in the context of classical doubles expanding in an ambi- 
ent medium with density profile p = p{r). These equations can be derived by (i) balancing the 
momentum flux of a relativistic jet by the ram-pressure of the ambient medium spread over some 
area A^, possibly larger than the jet cross-section Lj/c = pv'^A^, where Vh is the advance velocity 
of the jet head and Lj is the jet kinetic power; (ii) setting the lobe's sideways expansion velocity 
equal to the speed of the shock driven by the overpressured cocoon with internal pressure p in the 
surrounding medium, Vq = (p/p)^^^; and (iii) assuming that all the energy transported by a pair of 
jets during the source's lifetime t is transformed at the jet head (terminal shock) into the cocoon's 
internal pressure, pV = 2 (7 — 1) Lj t, where V is the volume of the cocoon and 7 = 4/3 is the 
adiabatic index of the ultrarelativistic cocoon's fluid. Introducing the source linear size LS and its 
transverse size Ic, one can therefore write 



Li = cp{LS)v^A]^ , p = p{l^)v^ , 3pV 
_ dLS _ die dV _ 

^ dt ^ ^ dt ^ dt 



2Ljt 



27r l1 



(1) 



For a given jet power Lj, source linear size LS", and ambient medium density profile p{r) hereafter 
assumed to possess King-type form, one also has to introduc e some additional scaling between the 
model parameters. Here we follow iKawakatu &: Kind ( 2000 ) with l"^ oc and fix 5 = 1 in order 
to reproduce the initial (ballistic, or 'ID') pha se of the jet propag ation into a uniform ambient 
medium as found in the numerical simulations of I S check et al.l (|2002l ). We also restrict our analysis 
to young GPS sources, which evolve in the central plateau of the galactic gaseous halo and thus have 
LS < Tn ~ 1 kpc. Therefore , we set the ambient density profile as p = mp no with no ~ 0.1 cm~^ (see 
Mathews &: Brighentil l2003l ) . Such a choice gives vy^ oc LS^, Vc oc LS"^/"^, p oc LS^^, Ic oc LS^/"^, 
V oc t oc LS, and A^^ oc LS^. Note, that the source linear size LS is defined here as the jet 
length, i.e., as the separation of the terminal hotspot from the core. 

In order to fix some oth er model parameters, we recall here several observational findings. 
Katz-Stone &: Rudnickl (jl997l ) analyzed the radio emission of the two CSS sources 3C 67 and 
3C 190 and found that the spectra of their lobes are consistent with relatively young source ages of 
t ~ 10^—10^ yr for the equipartition magnetic fields of i3 ~ 1 mG. This implies that the hotspots' ad- 
vance velocities, ?;h ~ 0.3c, are significantly (an order of magnitude) higher than the analogous val- 
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ues found for classical doubles. iMurgia et al.l (|l999l ) performed spectral ageing studies for a number 
of other CSS sources and found that, assuming again energy equipartition, the resulting source ages 
are indeed < 10^ yr and imply the average advance velo cities Vh ^ 0.3c. Such high velocities were 
i n fact detected direct ly, first in GPS sources 0710+439 (jOwsianik ConwavlEggsl ^ and 0108+388 



([Owsianik et al.l Il998l ). Since then, several radio obse rvations of hotspots in many GPS/CSO 



objects have confirmed repeatedly that ^ 0.2 h (jTavlor et al 



Giroletti et al.ll2003 



Polatidis fc ConwaylboOS : Gugliucci et al 



2005 



200C: 



200C- 



Tschaeer et al 
Nagai et al .I bood :! iGugliucci et al 



200?! : Iluo et ahlboO?! ). This value should, in fact, be divided by a factor of 2, because typically 
in the literature the relative hot-spot separation is discussed, i.e., the separation of one hot-spot 
relative to the other, and not the separation of the hotspot relative to the very faint (and therefore 
often undetected) core. In general, the established agreement between kinematic and spectral ages 
for the GPS and CSS populations (< 10^ yr for GPS and ~ 10^ — 10^ yr for CSS objects) supports 
approximately the fulfillment of energy equipartition within their lobes, at least in the majority 
of studied sources (magnetic fields ~ 10 mG and ~ 1 mG for GPS and CSS classes, respectively). 
In addition, it implies — by means of ram-pressure arguments — an ambient gaseous/interstellar 
medium (ISM) density tt-ism ^ 1 cm~^ on scales between a few and a few hundred pc (in agreement 
with the value anticipated here). 

Constant advance velocity of GPS/CSS objects in the simple model presented here, oc 
LS^, is in good agreement with observations, and thus, we fix hereafter vi^ ~ 0.1c. Another 
constraint is provided by the established (approximate) minimum power condition. In general, 
one can parameterize the magnetic field energy density in the expanding lobes as Ub = t/b p, with 
r/B < 3. Hence, with the model pressure 



P 



Lj mpUQ Vh 



1/2 



^ 5.3 X 10"^ L^/fg L^iVo erg cm 



1-1 



(2) 



where LSiqq = L5'/100pc and -Lj,45 
is expected to scale like 



Lj/10 ergs , the magnetic field intensity within the lobes 



1/2 .1/4 



-V2, 



(3) 



This is consistent with the equipartition values B ^ 1 — 10 mG typically obtained for the CSS 
(LS" ~ 1 - lOkpc) and GPS {LS ~ 0.01 - 1 kpc) objects, if r]]l'^L]^^t^ > 1. Note the comfortably 
weak dependence of the model magnetic field B on the jet kinetic power Lj. Note also that the 
agreement between the ages of GPS sources derived by means of the spectral ageing analysis and 
the dynamical one is in fact expected in the presented model since, in its framework, the age of the 
source during the GPS phase of the evolution is simply 



t 



-1 



LS w 3.3 X lO'^ LS 



The other model parameters can be evaluated as 



Ah 



100 yrs . 



2.2 X 10^°Lj,45cm^ 



(4) 



(5) 



-7- 



V 



3 TTipno 



^^^1 ^ 3.7 X lO^o LS\'l cm , 

1/2 



LS^ 1.3 X 10^2 ^i/J, ^52^^ cj^3 



Gvr rUpUQ J 



\ 1/4 

] L5-V2«l.8xl09Ly,U5 



-1/2 

'j,45 -^^100 



cms 



(6) 
(7) 
(8) 



Interestingly, the sideways expansion velocity Vc, again weakly dependent on the jet power, 
is rather high but consistent with the outflow velocities of the line-emitting gas (believed to be 
pushed out and accelerated by the expanding lobes), fout cms , observed at optical fre- 

quencies (with lumiiiosities of about ~ 1 0'^^ — lO'^^ ergs~^ for individual lines) in many CSS objects 
(jPe Vries et al.lll999l : lO'Dea et al.l 12002 ). In addition to this, gaseous outflows in host galaxies of 
GPS/CSS sources may be manifested as blue-shifte d absorption-line s ystems at UV frequencies. 
This, in fact, was observed in the CSS quasar 3C48 (jGupta et al.ll2005l ). indicating outflow veloci- 



ties of, again, Vont ~ 10^ cm s , driven by interaction of the expanding jets/lobes with the kpc-scale 
gaseous environment. We note that, since the Vc evaluated here may be higher than the expected 
sound speed in the external medium, Cg = (5A;T/3mp)^/^ ~ 3.7 x 10''cms~^ (for the anticipated 
temperature T ~ 10^ K characterizing the hot phase of the gaseous environment), a bow shock 
may be expected to form around radio lobes of GPS sources, with a possibly high Mach number 
A^sh = (3p/5noA:r)i/2 ~ 48 4/4^5 L^-^' . 



2.3. Synchrotron Emission 

2.3.1. Synchrotron Luminosity Evolution 

With the electron energy distribution injected from the terminal jet shock to the expanding 
lobe and modifid thereby by the adiabatic and radiative cooling effects, iVe(7), one can express 
the lobes' synchrotron luminosity as Lgyn = (4c(TT/3meC^) 7(7) U-q Ue V, where 7(7) = (7^)/(7) = 
J 'y'^ Nc{'j) dj / J J N(.{'y) dj and the electron energy density is simply Ue = rUeC^ J 'yNe{'y)d'y. 
Assuming further that Ue = rj^p with r/g ^ 3 (i.e., that the jet electrons shocked at the terminal 
hotspot, possibly in rough equipartition with the magnetic field and relativistic protons, provide the 
bulk of the lobes' pressure) and that the lobes' electron population does not change significantly its 
spectral shape during the GPS phase of the expansion (see § 3.1), the synchrotron luminosity turns 
out to be constant with time and independent of the source linear size, Lgyn (x UbUsV (x LS^, in 
particular 

1 /2 

^^-^ = ^Br?e/(7)4^'- 

« 1.2xl0^^r]Br]efi-f)L^%eigs-\ (9) 
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Such an evolution is expected to hold only if the source is not older than 10^ yr (i.e., the jet is 
in the initial, ballistic evolution phase) and if the ambient medium density profile can be approx- 
imated as being constant (i.e., if LS < Ikpc, the typical core radius of the gaseous medium in 
giant ellipticals). However, for t > 10^ yr, the above scaling breaks down, and the synchrotron 
lur ninosity is expe c ted to decrease w i th iii cr easing linear size of the source (see the discussion 



m 



BegelmanI Il996l : iKawakatu &: Kind l2006l ) . Note that equation 9 implies the interesting con- 
straint 1 < ??B ^^Ja5 f{l) ^ 10^ the presented model because the observed radio luminosi- 
ties of GPS sources (generally selected so far at intermediate/high flux density) are in the range 
Lsjn ~ lO''^ - lO^'S ergs-i. 

The scaling of Lgyn given above corresponds to a constant electron injection provided by the 
terminal hotspot, no absorption of (radio) synchrotron photons within the lobes, and to the fixed 
spectral function 7(7) during the source's evolution (see in this context §3.1). Note also that 
the monochromatic synchrotron power produced by electrons with a given Lorentz factor 7*, for 
which absorption effects can be again neglected, goes like [vLy\j^^^i,2 cc U^UeV cc LS^. On the 
contrary, the monochromatic synchrotron power measured at a fixed observed frequency i'* (and 
thus produced by the electrons with different energies at different evolutionary stages, due to a 
change in the lobes' magnetic field intensity) scales as [u^L^,*] oc i?('^~3)/2 (JeV oc LS^^"*)/^ for 
a power-law electron energy Ng^'j) cx 7"^ since the monochromatic luminosity can be written as 

WL^hyn = \caTVUB [7^A^e(7)] 



2caT VUbU^ /47rmeCi/^ (^-"^/^ 



3 nieC^ J ^ d'-f \ 3eB 



for iVe(7) oc 7-" . (10) 



2.3.2. Absorption Effects 

The observed turnover in the radio spectra of GPS/CSS sources is their main characteristic. 
It was proposed that it is due to either synchrotron self-absorption (SSA) or free-free absorption 
(FFA) by an inhomogeneous screen of dense ambient matter. The observed spectral indices below 
the peak frequency are usually aiow > —2; in some cases, they are close to the standard value 
—5/2 predicted by the homogeneous SSA model while in other cases, they are even consistent 
with the exponential cutoff predicted by the simplest version of the FFA model. The variety of 
the low-frequency spectral indices thus indicates inhomogeneity of the absorbing r nedium and/or 



superp osition of several emission components with different physical parameters. iDe Vries et al 



(|l997l ) found that the average spectral indices for the analyzed sample of GPS/CSS sources are 
oiow = —0.51 (±0.03) and ahigh = +0.73(ib0.06) below and above the peak frequency, respectively 
and that, in addition, the values of Ohigh are characterized by a very broad distribution between -1-0.5 
and +1.2. They also claimed a flat spectral plateau between fp and 2 x i/p in the template GPS/CSS 
spectrum, with average power-law slope apeak = +0.36(±0.05), ascribed to the broadening of the 
spectral peaks by the sources' substructure. 
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An implication of the evolutionary models for GPS/CSS sources with strong negative lumi- 
nosity evolution (i.e., t he models assum ing steep ambient medium density profile pismif) ^ 
with P > 1, like that of iBegelmanI 1 19961 ) is that the young sources evolve on the fp — LS plane as 
Up oc LS~^ with X > 1 if the spectral turnover is due to the SSA process. Thus, the GPS/CSS 
sources do not evolve along the observed Up oc ^g-O-S S ^^^g 'leave' the Up — LS plane when 



their radio powers decrease enough with increasing LS (lO'Dea &: Baumlll997l). This could possibly 



explain the observed scatter in the fp — LS correlation. On the other hand, iBicknell et al.l ()1997l ) 
successfully reproduced the observed Up oc LS~^'^^ dependence in the framework of a model in 
which the spectral turnover is due to FFA by a clumpy /filamentary multi-phase ISM, modified 
(ionized) by the p assage of a bow shock due to expanding radio lobes (see also in this context 
Kuncic et al.lll998l ). However, this model requires very special parameters for the ambient medium 
(its high density, in particular) not consistent w ith the most recent observations for the majority 



of sources (see the discussion in iBegelman 



19991). Th e othe r promising possibility left is therefore 



the 'engulfed cloud' scenario proposed bv IBegelmanI (jl999l ). in which the neutral clouds of ISM 
penetrating the expanding radio lobe and photoionized by the nuclear radiation are responsible for 
the spectral turnover of GPS sources due to free-free absorption of the radio photons. 

Can the SSA effects be responsible for the observed spectral turnover of young radio sources 
in the framework of the dynamical model proposed here? To investigate this, we note that the 
characteristic SSA frequency can be found from the equation Tjf^ = 1, where t^^^ = LS is the 
optical depth for the synchrotron self-absorption process described by the absorption coefficient 
^ssa -^^^j^in the uniform medium of the spatial scale LS. Since, in the case of a power-law electron 
energy distribution, this coefficient is n'^^ oc A'^g B^'^"*"^-'/^ 1/"^*"*"^^/^, the characteristic (peak) SSA 
frequency is fgsa oc LS~^, with x = (s-|-2)/(2s-|-8) = 0.3 — 0.36 for s = 1 — 3. This is flatter than the 
observed distribution of the turnover frequency in GPS/CSS sources, Vp oc LS~^'^^. On the other 
hand, as noted before, the observed Up — LS distribution may be shaped by some additional factors 
related to the sources' luminosity evolution and to the incompleteness of the samples considered 
and thus may not reflect directly the evolution of the peak frequency in a single source. Note, in 
this context, that in the framework of the discussed model (see equation 10) the synchrotron power 
at such an SSA frequency scales as [z^ssa-^^vssa] o<i B^^~^^/^ Ussa. ^''^^ UbUsV oc LS^, with y = 0.2 — 
for s = 1 — 3; i.e., it is expected either to increase slightly or even to remain constant with increasing 
linear size of the source. Thus, it seems unlikely that the evolving young s ources 'leave' quickly t he 



flux-limited fp — LS plane due to a decrease in their peak luminosity (cf. iQ'Dea &: Baumlll997l ). 



In addition, SSA effects are expected to manifest themselves at relatively low frequencies. To 
illustrate this, let us assume that the initial electron energy spectrum injected from the hotspots 
to the expanding lobes is, on average, of a power-law form with the 'standard' spectral index s = 2 
(although the obtained results hold approximately for a broader range of 1 < s < 3). In such a 
case, one gets Ue = rrieC^ Nq ln[7max/7min] « 10 mgC^ A^o, and 



3e^ 



4tt mic 



3„3 



No y- 



(11) 
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Hence, the condition r^*"^ 



1 gives the critical SSA frequency 
'0.148 nip np vi, rjc tje Lj 



Wnm^c^LS 
l/3„l/3 .1/3 



1/3 



0.3 r/,^/^ 77b 



1/3 



GHz. 



(12) 



Because r/B, % ^ 3 and the expected jet kinetic power is Lj < 10^^ ergs~^ (see below), this frequency 
is much lower than the turnover frequency observed in GPS/CSS sources, allowing us to neglect 
SSA effects in the following discussion. 

The above conclusion is in agreement with several observational supports for t he free-free 



absorp t ion process shaping the in verted spect r a of GPS/CSS objects. For example, iPeck et al 



( 


2000), and 


Marr et al. 


( 


2001) 



OQ 208 and 0108+388, respectively, based on investigations of spectral index maps and hence on 
differences between low-frequency spectra in different parts of the sources. Their modeling implies a 
non-uniform absorbing gaseous medium with high average number density (~ lO^cm"^), ascribed 
by the authors to, e.g., a clumpy torus-l ike structure o f the o bscuring material with extension 
< 100 pc. Supporting this interpretation, iKameno et al.l (j2003l ) showed that, on average, type 1 
(Seyfert 1, quasars) and type 2 (Seyfert 2, radio galaxies) GPS sources have different characteristics 
of the absorption features (asymmetric and symmetric distributions with respect to their nuclei, 
respectively) , con sistent with the FFA process associated with anisotropically distributed obscuring 
material. Finally. iMutoh et al.l (|2002l ) noted that the SSA model implies a change of the polarization 
angle by 90° across the spectral peak, i.e., between the optically-thick and o ptically-thin parts o f 
the continuum. Since no such changes were observed in the sample considered, iMutoh et al.l (|2002l ) 
argued for FFA effects playing the major role. One should be aware, however, that the GPS quasars 
may not be in reality compact/young sources, as emphasized in § 1 of this paper, and therefore 
that they may be intrinsically different from the truly young GPS radio galaxies. 



Let us therefore discuss in more detail the absorption model proposed bv iBegelmanI (jl999l ). In 
this model, dense hydrogen clouds of ISM present at pc— kpc distances from the centers of young 
radio sources and engulfed by their expanding lobes are photoionized by the active nuclei, causing 
free- free absorption of the lobe radio emission. Such clouds may be naturally identified with the ones 
producing narrow-line emission and also Hi absorption lines. Very broad absorption lines are indeed 
often detected in GPS/CSS objects (at much higher rates than in extended radio g alaxies) with neu- 



tral hydrogen column densi t ies from iVHi 5-, 10^ ^ cm ^ down to A'^hi ^ 10^^ cm ^ ( Vermeulen et al 



20031 : iPihlstrom et al.l 120031 : iGupta et al.l 120061 ). T he evaluated column densities were claim ed to 
anticorrelate with the source sizes, A^Hi oc LS-^-^^ (IPihlstr5m et al.ll2003l : iGupta et al.ll200(il ). Pre- 
viously, it was speculated that the appropriate absorbing medi um possesses a, torus-like distribution 



However, detailed studies of a few objects performed to date (jLabiano et al.ll2006l : IVermeulen et al 



20061 ) indicate that the Hi absorption lines are associated with optical emission lines and thus 
arise most likely in the atomic cores of NLR clouds interacting with the expanding radio source. 
Interestingly, the intensity ratios of lines produced by spatially resolved NLR clouds in nearby 
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Seyfert galaxies imply a densi ty decrease with d i stance from the ion i zing source (galactic center) 



PC r "•, with 1 < n < 2 (see iKaiser et al 



Kraemer et al. 2000: Mundell et al 



2003 



200C 



Nelson et alj l2000l : iKraemer fc Crenshaw! l2000l : 



Bradley et al.l 12004 ) . This would then be consistent with 



the noted anticorrelation of A^hi with LS if Hi absorption is indeed due to NLR clouds and not to 
the hot phase of the ISM gas (which then r nay obey a King- type distribution with central plateau 
"-ISM oc const, as assumed in this paper; cf. iPihlstrom et alJl2003l ). 



It is also interesting to note that GPS sources exh ibit in general (ve ry) low polarization of 

, while radio continua 
Such low polarization 



their radio fluxes (on the level of a few percent, if any: iQ'Dea et al 



of CSS objects are polarized a bit more strongly (jStanghellini et al 



1991 



2001 



is most probably due to Faraday effects. The rotation measure observed has a v ery broad scat- 
ter in the GPS/CSS sample, from very large, RM > lO'^radm-^ (|NanetalJ l2000|) , to very small, 
RM < lO^radm-^. In some cases, RM is very different for the lobe and the counterlo be in a single 



source, which indicates large asymmetries in the gaseous environment on kpc-scales (jJunor et al 



19991 ). Faraday scree ns seem to be assoc i ated predominari t ly with the optical line-emitting clouds 



interacting with jets (jCotton et al.ll2003l : iFanti et al 



2004 



Cotton et al.l 12006 ). Indeed, the jets in 



CSS quasars e xhibit more complex an d distorted radio morphologies than the jets in flat-spectrum 
radio quasars ([Mantovani et al.lll998l ). suggesting significant jet- ambient medium interactions tak- 
ing place somewhere within the NLR. Also, some other morphological and polarization properties 
often suggest an asymmetric distribution of the ambient medium surround ing CSS objects, p o ssibly 



resulting in an intrins ically asymmetric structure of their lobes and jets (|Saikia et al. 
Saikia fc GuptalboO^I ^. 



1995 



2001 



In a scenario where the engulfed photoionized NLR clouds are responsible for the Hi absorption 
lines and the d istortion of the radio structures but are unable to confine or frustrate the jets 
(jBegelmanl 1 1 999l ) . the dominance of FFA effects by those clouds in generating the spectral turnover 
in the radio spectra of GPS and CSS objects would imply a phenomenological relation between the 
peak frequency and the source's linear size, Up ~ 2.7 LS^qq^^ GHz. The synchrotron luminosity at 
the peak (turnover) frequency is thus expected to scale with the size of the source as [i^pL^^] oc 
^(s-3)/2 j,j3~^)/2 UbUcV (X LSy, with y = (-0.15) - for s = 1 - 3 (equation 10); i.e., it should 
decrease only very slightly with increasing LS or even remain constant during the GPS evolution 
phase. Note that the low-frequency radio continua are then expected to be of the form L^^^^ oc 
^^/^ because, in the describe d engulfed-cloud model, the appropriate absorption coefficient 



,2-(s 



scales as oc ( Begelmanlll999l ). Therefore, the observed scatter in radio spectral indices below 
the spectral turnover may result not only from the non-uniformity of the lobes and NLR but also 
from the internal scatter in the low-energy electron spectral index s in the lobes of GPS objects. 
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3. Broad-Band Spectra 

3.1. Electron Energy Distribution 

It is typically assumed that the electron energy distribution formed at the jet terminal shock 
and injected to the lobes is of a simple power-law form Qij) oc 7^*. Such an injected spectrum 
undergoes further radiative and adiabatic energy losses within the expanding cocoon. For the GPS- 
phase of the evolution described in § 2.2, the lobes' expansion is given by V/Vq = -kII/A]^ oc LS, 
and the cooling (magnetic and radiative fields) changes with the source size LS. Assuming that 
synchrotron emission dominates the radiative losses (as is the case for all but the most compact 
sources, as discussed in the next section § 3.2), the resulting electron energy distribution for a given 
source linear size LS = v^t with Ub oc LS~^ can be found from equation A8 of Appendix A, which 
reads 

A4(7) = y^^Q[7(i^5')] j{LS'fLS"/'dLS', (13) 

where Meij) = N^ij) V. In the integral above, 

1 + C3 70 [LS-^^^ - LS'-^/^) 

(see equation A6), where the constants C2 ~ 2.5 x 10^ -^j^ 45 cm^/^ and C3 1.34 x 10^ t/b -^j^ 45 cm^/^ 
follow from the equations 2, 5, and 6. For a given power-law source function Q{'y) = Kq^~^ with 
s > 2, one can obtain further 



s{s + l){s- 1) (s + 2) 

s(s2-l) 3s(s-l) 6(s 



A4(7) 

s ys -\- L) \s — J 

rs(s2-i) 

""1 {u + l) + {u + lf +(a; + l)3 + (u; + l)4/' ^^^^ 
where, after integrating, 70 was replaced back with jLS^^^/ [c2 + C3 7 (l — {LS / LSq)^^^)] , and 

C2/C3 _ 7cr ^ „„n -1 r-V2 -1 

= V 200 Lj 4^ 7 . (16) 



CO 



7 7 



Equation 15 implies that the electron energy distribution in the lobes of GPS sources, in the 
case of a single power-law injection oc 7~'' with s > 2, is expected to be of a broken power-law 
form, 

^ . . _ AOt) _ ^^e^ / 7"' for 7<7cr 

"^^^^^ - ~ ^ M 7. 7-^-^ for 7 > 7c: • ^''^ 



Note that the critical break 7cr does not depend on the linear size LS and hence neither on the 
function 7(7) = / 7^ ^eil) dl/ J 7 A^e(7) d'y ~ 7cr*/(3 — s). In addition, since V oc LS^ (equation 
7), the normalization of the electron energy distribution for a given LS scales as LS~^. Similarly, the 
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electron energy density Ue = nieC^ J 7 A'e(7) oc LS~^, ensuring that the ratio Uc/Ub is constant 
during the GPS phase of the lobes' evolution. These scaling relations validate the discussion 
presented in § 2.3 and allow one to fix the normalization of the electron injection function, Kg, 
through the assumed relation Ue = ijep. One may finally find that for a single power-law injection 
with s ~ 2 — 3 considered in this paragraph, the total (unabsorbed) synchrotron luminosity of GPS 
radio galaxies is expected to be roughly 0.1 Lj < Lsyn ^ Lj (see equation 9). 

In the discussion above we have assumed the injection of a single power-law electron energy 
distribution within the terminal hotspots of GPS sources. However, sorae rece nt studies indicate 
that the situation may be more complicated. For example, IStawarz et al.l (j2007l ) found that in the 
case of the archetype FR II radio galaxy Cygnus A, the electron spectrum can be approximated as 
a broken power-law Qi'j) oc 7"*^ for 7 < 7int, and Q{'y) oc 7"*^ for 7int < 7, with 7int ~ mp/m-e, 
si ~ 1.5 and S2 > 3. Such a form is, in fact, expected in the case of cold protons carrying the 
bulk of the jets' energy in powerful radio sources (at least on > pc scales) due to the nature of 
the particle acceleration process taking place at t he mildly relatiyistic terminal shocks dynamically 
dominated by the protons (see the discussion in IStawarz et al.l 120071 ) . If this is also the case for 
young radio sources, then one can expect the electron energy distribution within the lobes of young 
radio galaxies, which can be evaluated as 



iVe(7) = 7- 



LS 



Q 



4/LS 



uj 7 



(lj + l)a;V3 _ 1 



oj 7 



fa; + l)xV3 _ 1 



dx 



(18) 



(see equation 13), to be peaked on the 'y'^N^'y) — 7 plane around the electron energy ■jcr, as long 
as si < 2. Hence, one gets 7(7) ~ 7cr, and therefore, the expected total synchrotron luminosity 
Lsyn ^ 0.1 rye ij, similarly to the case of a single power-law injection discussed above. Thus, 
the proposed model with the anticipated r/e ^ 1 is in agreement with the observed values of 
10^^ — 10^^ ergs^^ for different injection conditions and a jet luminosity in the range Lj ~ 
lO^'^ergs"^. This range is indeed as expected if GPS sources are progenitors of FR I and 
FR II radio galaxies and implies that > 10 % of the jet kineti c power is dissipated for the synchrotron 
emission of the GPS lobes (as assumed bv lDe Young||l993l ). 



Lsyn ~ 



10' 



,43 



The evolution of the electron energy distribution within the lobes of GPS sources, as given 
by equation 18, is shown in Figure 1 for different jet powers (Lj = 10^^, 10^^, 10'^^, and lO'^^ erg 
s~^), and different source linear sizes {LS = 100 pc and Ikpc; thick/upper and thin/lower lines, 
respectively). In the figure, two different source functions are considered for illustration, namely 
single power-law (3(7) oc 7"^-^ (dotted lines), or broken power-law Q{'y) oc 7"^-^ for 7 < 7int and 
(5(7) oc 7^^ for 7 > 7int (solid lines). In both cases the minimum and maximum electron Lorentz 
factors are 7min = 1 and 7max = 10^, while the normalization of the injection function is evaluated 
through the condition Ue = r]eP with rje = 3 and t/b = 0.3. Vertical dotted and dashed lines in 

1 —1 /2 

the figure indicate critical electron energies 7cr = 200 7?B'Lj 4^ and Tint = mp/nie, respectively. 
As shown, the normalization of the electron energy distribution decreases with increasing source 
linear size LS roughly as LS^^ while the spectral continuum steepens at 7 > 7cr when compared 
to the injected one. Note also that, in the case of high jet luminosities and small source linear 
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sizes, the peak in the electron energy distribution for a broken power-law injection may be slightly 
higher than 7cr- Thus, the total synchrotron luminosity may exceed the value 0.1 rjcL^ for the most 
powerful and compact GPS sources. 

A more complex form of the electron source function Qij) introduces several other interesting 
features in the model. In particular, with an intrinsically broken electron spectrum, further modified 
by the adiabatic and radiative energy losses, one should expect to observe a curved or multiply bro- 
ken synchrotron continuum. That is because the synchrotron spectrum in such a case is shaped by 
three different critical frequencies: one related to absorption effects, fp ~ 2.7 LS~^'^^ GHz, and the 
other two related to the critical electron energies discussed above, namely i^cr = 3 e-B 7^j./47r rrieC ~ 
0.54 ?7b^^^ LS^qJ'^ GHz, and Umt = 3 eB jf^^^/Air rricC « 70 r]^"^ -^j 45 -^"^loo^^ GHz. All these are 

interestingly clustered around GHz frequencies; some spectral plateau around this range is thus to 
be expected. This kind of a spectral feature is indee d observed in many p articular objects and is 



also present in a template GPS spectrum derived by lDe Vries et al 



lany p a 



3.2. Photon Fields Within Lobes 



At near-UV frequencies, GPS/CSS sources, like classical doubles, exhibit complex spectra 
comp osed of nebular continuum, nuclear light (bo th direct and scattered), and a starburst compo- 



nent (Tadhunter et al 



2002 



De Vries et al.l 120071) . An additio nal non-thermal con tribution from 



j ets a nd compact lobes was also anticipated (jBegelmanI Il999l ). As suggested bv iLabiano et al 



(1200a), the observed extended UV luminosities of several GPS/CSS objects, being in the range 
10^^ erg s~^ , are correlated neither with the sources' radio powers nor with their linear 



10 



sizes. 



40 



Labiano et al.l (120081 ) argued that this emission, produced predominantly by relatively young 



stars, results from bursts of star formation that took place at the time of or before the formation 
of compact radio structures. Moreover, the detected (in some cases) UV component aligned with 
the radio axis may also indicate star formation enhanced or triggered by the expanding radio lobes. 
Here we concentrate on the UV photons provided directly by the active center, assuming that nu- 
clei of GPS objects are intrinsically similar to the ones observed in quasars and Seyfert galaxies. 
Therefore, we assume that the bulk of the radiative output of the optically thick accretion disk 
is emitted at UV frequencies (thus forming the characteristic 'big blue bump') very close to the 
central engine, with iii trinsic luminosities of the order of Luv ~ (10^^ — 10^^)ergs~^ (see, e.g.. 



Koratkar &: Blaeslll999l ) and thus dominating over the other, extended UV ph oton fields mentioned 



abov e. Such strong emission is in fact observed directly in many GPS quasars (ISiemiginowska et al 
20051 . also in prep.). In the case of GPS radio galaxies, however, the UV disk emission toward the 
line of sight is likely to be absorbed by obscuring dusty tori. 

Note that, in the framework of a model ascribing the spectral turnover of GPS sources to 
free-free absorption by the engulfe d NLR clouds, it is the UV disk emission which is required to 
photoionize the absorbing matter (jBegelmanI Il999l ). In addition to this, the UV disk radiation 
provides an important source of seed photons with energies hi'Q ~ 10 eV for the inverse-Compton 
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(IC) emission of ultrarelativistic lobe electrons. The volume- averaged energy density of this photon 
field, as a function of the distance from the central engine (and hence the source linear size LS), 
can be simply estimated as 



1 

V 



L 



uv 



4ttLS^c 



dV 



L 



uv 



10 ^ L\]y^iQL S^QQ dig c-m 



2ti LS'^c 

3 



In 



LSn 



(19) 



3 and Luv, 46 = Lyy /W^^ ergs ^. This is comparable to or 



where we put In (L5'max/i5'min) 
slightly lower than the magnetic field energy density within the lobes, Ub 
(see equation 2), for the expected values -^^j,45 > 0.1 and Luv,46 < 10. 



m w-^lI%ls- 



100 erg cm ^ 



Some part of the strong direct disk emission in powerful AGNs — typically more than 10% — 
is expected to be reprocessed by the obscuring matter (dusty tori) and re-emitted at FIR-to-NIR 
frequencies. Several observations confirm th e presence of such a s pectral feature in the radiative 
outputs of young radio galaxies. For example. iHeckman et al.l (jl994l ) showed that GPS/CSS sources 
have the same MFIR strengths as extended sources with comparable radio powers and redshifts, 
i.e., (L.ejQ urn) ^ 3 X 10^^ e rg s^-^ for (^sgHz) ~ lO^'^ergs"^ in the case of GPS/CSS radio galaxies 
(see also iHes et al.l Il995l ). We note, however, that the extended sources pr obably have lowe r 



efficiency in converting the jet kinetic power to synchrotron radio emission. iFanti et al.l (l200d ) 
confirmed that the dust emission dominating the radiative output at FIR frequencies is similar in 
the GPS/CSS population and regular radio galaxies, both in luminosity and temperature. This 
implies that the young radio sources do not contain more (or less) dust than the extended ones: on 
average, the FIR luminosities of these objects are (^^fir) ~ 3 x 10^"^ ergs~^ and can be modeled by 
a two-tempe rature dus t distr ibution with a total mass of ~ 10^ Mq. These findings were recently 
confirmed bv lShi et al.l (|2005l ). who emphasized similarities in MFIR emission between CSS objects 
and regular radio- loud quasars or powerful radio galaxies. 

The energy density of the photon field due to the dusty torus at a dist ance r from the galactic 
center can be estimated as (see lSikora et al.l 120021 : iBlazejowski et al.ll2004l ) 



Umir) 



1 



Anrjc 1 + {r/rdY 



(20) 



where Ljr ~ Luv is the expected torus luminosity, and = {L\jy /An asBT^)^^"^ is the charac- 
teristic (minimum) distance of the circumnuclear dust with temperature ~ 10'^ K. Since is 
supposed to be much smaller than the spatial scales considered here, r 
restrict the analysis to LS > r^, thus obtaining the volume-averaged 



■^uv 46 P^' ^"^^ ^^"^ 



Um ~ 10" 



-^-IR, 46 -^-'S'loo erg cm 



-3 



(21) 



with LiR,46 = LiR./lO'^^ergs-^ The mean value of the IR photon energy is hereafter assumed as 
~ 10"*^^ Hz. That follows from the fact that the dust temperature is oc r^^/^, while the dust 
density is inversely proportional to some high power of r. Therefore, the IC emission taking place 
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at larger (r > rj) distances, as considered here, is mainly related to the hotter dust, i.e., to the 
MIR frequencies of the target photons. 

The photon fields f/uv f^iR evaluated above scale linearly with the nuclear luminosity 
Luv- So far, we have assumed a relatively high value of Luv = 10^^ ergs~^, as is appropriate for 
the nuclear spectra of quasar sources. It is not clear, however, whether the nuclei of all GPS/CSS 
radio galaxies are indeed similar to the quasar ones. The emerging consensus is that powerful radio 
galaxies of FR II type — those possessing high-ionization broad and/or narrow lines — are in fact 
misaligned quasars, with the UV accretion-related emission obscured toward the line of sight by 
dusty tori (|Barthellll989l ). This is supported by the fact that the obscured nuclear UV component is 
often observed in these objects indirectly via polarized scattered light or via intense re-emission of 
the obscuring matter at MFIR. In the framework of such a unification scheme, radio-loud quasars 
are eventually brighter in IR when compar e d to powerful FR lis only due t o the contribution o f 



the beamed jet emission (IHaas et al 



2004 



Shi et al. 



2005 



Ogle et all |2006|; IClearv et all |2007|). 



Thus, powerful GPS/CSS sources which are precursors of powerful FR II radio galaxies are not 
expected to lack the strong accretion-related UV and IR emissions discussed above. The situation 
is, however, less obvious in the case of young radio sources that are precursors of FR I radio galaxies 
or of weaker FR II objects with only low-ionization emission lines. 

The radiative output of active nuclei in FR I radio galaxi es at optical frequenc ies is most 
likely dominated by the synchrotron emission of unresolved jets (IChiaberge et al.ll2000l). They are 
generally believed to lack strong accretion-related optical/UV radiation ( but see iMaoz 2007), i n 
agreement with the weak nuclear obscuration advocated for these sources (IChiaberge et al.ll2002l ). 
Nuclei of lower-power FR II sources are very often similar to nuclei of FR I radio galaxies, with 
only low-ionization emission lines present, optical flux most likel y dorn i nated by the jet emission, 



and relatively moderate MFIR emission (jChiaberge et al. 



200C 



2002 



Ode et al.ll2006l ). Thus, 



nuclear obscurat i on in these objects is believed to be only moderate (although still present; see 



Haas et al. 



2005 



Hardcastle et al.ll2006l ): in other words, accretion-related emission is believed to 



be lower than in powerful FR Ils/quasars. Clearly, in the case of GPS/CSS radio galaxies that 
are precursors of these objects, kpc-scale photon fields might be lower than our model parameters 
Luv ~ -^^IR ~ 10^^ ergs~^. Even in this case, however, the synchrotron emission of the radio lobes 
provides a relatively intense photon field for the IC process. Its energy density can be evaluated as 



j^j- _ J syn dv 



(22) 



where the synchrotron luminosity (equation 10) is integrated over the whole frequency range with 
absorption effects included. Taking approximately / Ly^^y^dv ~ L^y^ as given by equation 9 (i.e.. 



ignoring absorption of the radio photons), as well as /(7) ~ 7cr (see section § 3.1), one obtains 



5.2 X 10-9 r?e L5ioo^' erg cm'^ 



(23) 



which is again less than or, at most, comparable to the magnetic field energy density U-q. 
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Yet another source of seed photons for IC scattering within the lobes of young radio objects is 
provided by the optical light of host galaxies. In general, the host galaxies of GPS / CSS sources are 
very simila r to the hosts of powerful (SCR) classical doubles when observed in NIR (jPe Vries et al 



19981 . l200d ) , being evolved ellipticals with some morphological indications of relatively recent merger 



events. Very often, they a lso ex hibit kpc-scale opti c al eni ission aligned with the main axis of the 
radio source. lAxon et al.l (j200Cl ) and iLabiano et al.l (j2005l ) found that the aligned optical features 
consist of at least three distinct components, namely an unresolved nucleus, line-emission aligned 
with and of similar extent as the radio lobes, and a weak component extending beyond the radio 
lobes. The unresolved optical nuclei are present only in the sources with broad emission lines, 
in agreement with the nuclear obscuration/unification scheme. The kpc-scale line-emitting gas is 
most probably due to the interaction of gas clouds with the expanding jets and lobes. Finally, 
the diffuse component extending beyond the radio lobes is most likely photoionized by anisotropic 
nuclear emission. Hereafter, we consider conservatively strictly the starlight emission of the hosts, 
which corresponds to NIR frequencies of fo ~ 10^^ Hz. For a given IZ-band luminosity of the galaxy, 
Lv, which is known to cor relate well with the c ore radius of the stellar distribution, rg, namely 
Lv/10'^^ ergs~^f« rg/lkpc (jde Ruiter et al.l l2005l ) . the starlight energy density may be evaluated 
roughly as 



star 



3Ly 



8 X 10"^° Ly erg cm" 



(24) 



where -Ly, 45 = Ly/lO^^ ergs ^, and we assumed approximately LS < r^. 



3.3. Inverse-Compton Emission 

In the Thomson regime of the IC scattering, the IC luminosities related to the external photon 
fields can be simply evaluated as 

(25) 

Here ?7rad stands for the energy densities of different seed radiation sources (C/uV; f^iR; or C/gtar), 
which are approximated as being monochromatic with characteristic seed photon frequencies uq as 
given in § 3.2. In the case of SSC emission, the monochromatic approximation is not correct, and 
the appropriate luminosity can be evaluated as 

WL^Uc^ ^(TT-^ / ^ ^'^o,syn [l^ N^il)]^^^^^;;/^^ di^o ■ (26) 

The resulting broad-band spectra are shown in Figures 2-3 for Ly = lO^^ergs"^, different jet 
kinetic powers Lj (i.e., 10^^, 10^^, 10^^, and lO^^ergs"^), different source linear sizes LS (i.e., 
33 pc, 100 pc, and Ikpc), and fixed lobe parameters r/B = 0.3, r/e = 3. As illustrative values we 
also take Luv = -^IR = lO^^ergs"^ for Lj > lO^^ergs"^, while Luv = -^^IR = lO^^ergs^^ for 
Lj < lO^^ergs"^, accordingly to the discussion in section §3.2. Different figures correspond to the 
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different source functions Qij), namely, to the single power-law injection with spectral index s = 2.5 
(Figure 2) or to the broken-power-law with si = 1.5, S2 = 3 and the break energy 7int = mp/nie 
(Figure 3). The absorption, radiative cooling, and adiabatic cooling effects in the electron energy 
distribution A^e(7) are taken into account as described in §2.3 and §3.1. Finally, we assume 
minimum and maximum electron Lorentz factors 7min = 1 and 7max = 10^ > respectively. Solid 
lines in the figures correspond to the synchrotron component, dashed lines to the SSC emission, 
dotted lines to comptonization of the IR torus radiation, dot-dashed lines to comptonization of the 
starlight photon field, and short-dashed lines to comptonization of the UV disk emission. 

Although the model predictions regarding the high-energy IC radiation may differ from source 
to source due to the scatter in the model parameters, one can outline some of the main results of the 
analysis. First, the evaluated SSC, IC/IR and IC/UV components are together very pronounced 
at ~ keV photon energies. One should note, at this point, that the emission from accretion disks 
in AGNs extends as well to X-ray photon energies (being produced by Comptonization of thermal 
optical/UV disk radiation within the hot disk's corona) and that this emiss ion may be strong enough 
to dominate the inverse-Compton X-ray emission analyzed in this paper. iKoratkar &: BlaesI (jl999l ) 
showed that, in the case of quasar sources, the 1 keV disk luminosity is on average ~ 10^'^ times 
lower than the UV-bump luminosity. This, for Luv = Iff^^ergs"^, is still high enough to dominate 
the X-ray lobes' emission evaluated above in all but the most powerful and compact GPS quasars. 
In addition, the contribution from unresolved relativistic jets to the total X-ray output may be 
also non-negligible in these objects. Therefore, the keV inverse-Compton emission from kpc-scale 
lobes is expected to be pronounced in powerful GPS radio galaxies rather than GPS quasars since, 
in the former, the direct X-ray emission is expected to be absorbed by the obscuring tori and the 
jet contribution is expected to be Doppler-hidden. In a subsequent paper, we investigate in more 
detail whether the IC lobe emission can indeed account for most of the X-ray fluxes detected from 
GPS radio galaxies, with the intrinsic, i.e. absorption-corrected. X-ray luminosities typically in the 
range Lx ~ 10^^ — 10^^ erg s^^ and the fitted X-ray spectral indices < ax ^ 1-5. Interestingly, 
the lower-energy tails of the SSC, IC/IR and IC/UV spectral components may account for (or at 
least contribute significantly to) the extended UV emission det ected in a number of GPS/CSS radio 
galaxies, with observed luminosities Lx ~ lO^'' — lO^'^ ergs~^ ( Labiano et al.l 120081 ). 



We note, in this context, that the initial low detection rate of GPS/CSS quasars at hard X-ray 
frequencies suggest ed t hat these objects were either intrinsically X-ray weak or heavily obscured. 
O'Dea et al.l (|2000l ) and iGuainazzi et al.l (|2004l ) argued for the nuclear (disk) origin of the observed 
keV photons and significant internal obscuration, with the required absorbing column densities 
being as large as A^h ~ 10^^ cm~^ in some cases. Comp lex X-ray spectra cora plicate the interpre- 



tation of typically weakly detected spectral features (see IGuainazzi et al.ll2004l . for the case of Mkn 
668). The most recent X-ray observatio ns seem to confirm that GPS /CSS sou rces are obscured 



rather than intrinsically weak in X-rays (IGuainazzi et al 



2006 



Vink et al.l 120061 ). The implied ob- 



scuration is consistent with the presence of nuclear tori, with inferred absorbing column densities 
ranging from Ah ^ lO^'^ cm~^ up to Ah ^ 10^^ cm~^, but only in a few cases. Such obscuration is 
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similar to that observed in extended and powerful rad io galaxies, suggesting a very similar nuclear 
environment in GPS and FRII objects. As noted bv lVink et al.l (|2006l ). in all of the cases where 
the appropriate comparison could be performed, the evaluated A^h values exceed substantially (by 
a factor of > 10) the absorbing column densities A^hi implied by the detected Hi absorption lines. 
Such behavior may, in fact, be expected, taking into account the different locations of the absorbing 
media involved (nuclear tori vs. NLR clouds) and the different states of the absorbing matter (total 
vs. neutral fraction of the gas). 

Interestingly, in some CSS objects, clear excess (flat power-law) emission was noted at hard 
X-rays and as cribed to emission f rom the inner parts of the radio outflow. This is the case in 
quasar 3C 48 dWorrall et al.ll2004l'). Sevfert 1 object PKS 2004-447 (jCallo et ahlbood l. and radio 
galaxy 3C 303.1 (jO'Dea et al.l 120061 ). We note that the first two sources are especially bright in 
the infrared, most probably due to emission from the obscuring torus, and that this radiation can 
provide the dominant seed photon population for efficient X-ray inverse- Compto n emission in the 



inner (pc-scale) portion of relativistic jets in quasar sources (see in this context iGuainazzi et al 



2004 



Gallo et al.ll2006l ). Such a scenario was, in fact, proposed (and successfully applied) for th e 



particular case of GPS quasar PKS 1127-145 (jSiemiginowska et al.ll2002l : iBlazejowski et al.ll2004l ). 



Another interesting finding of our analysis is that the inverse-Comptonization of the UV disk 
emission seems to be strong enough to be detectable by GLAST for many GPS radio galaxiesEl 
Note that the estimates presented in this paper are rather conservative because we have assumed 
roughly equal energy stored in the lobe magnetic field and in ultrarelativistic electrons, with the 
magnetic field intensity being a factor of only a few below the equipartitiqn value (r?R = 0.3, 



3, in analogy to the lobes of classical doubles iKataoka &: Stawarz 



2005 



Croston et al 



2005 



Ve 

and references therein). Any larger deviations from the minimum power condition t/b ^ Vc would 
increase the expected inverse-Compton emission for a given synchrotron power. Thus, one can 
conclude that young precursors of powerful radio sources will most likely constitute a very numerous 
class of extragalactic 7-ray sources that will be detected by GLAST. It is interesting to mention, 
in this context, that in the most recent sampl e of radio-loud quasars extracted from a combination 
of the SDSS (DR3) and FIRST surveys by iDe Vries et al.l (120061 ) . the majority of objects are 
characterized by a compact radio morphology. As pointed out by the authors, beaming effects on 
a parent population of extended sources cannot account for the observed compact-source counts. 
Thus, many of these objects have to be truly compact, and a fraction of them must be young 
and powerful GPS sources. Again, in the case of such compact GPS quasars, the IG emission 
from the unresolved je t s may do minate the lobes' radiative output at GeV photon energies (see 
Blazeiowski et al.ll2004l : iBailbood '). 



In order to quantify prospects for the detection of GPS radio galaxies in 7-rays, let us estimate 
roughly the high-energy segment of the expected IG/UV emission, which dominates the lobes' 



^At lower 7-ray energies (< lOMeV), however, t he IC emission fro m lobes discussed here may be overwhelmed by 
the 'cocoon bremsstrahlung' radiation discussed bv lKino et al.1 l|2007l ). 



- 20 - 



radiative output at these photon energies. Assuming a single power-law injection with relatively 
steep electron energy index s = 2.5, corresponding to the 'radiatively cooled' emission spectral 
index a = i (s - 1) + 0.5 = 1.25, one may find from equations 17, 19, and 25 the monochromatic 
IC/UV luminosity 



-0.25 



1042erg/s V^b/ V 10^^ erg/s y \imY>c) VIO^^ erg/s^ V 1 GeV/ ^ ' 

or the appropriate IC/UV energy flux 



10-12 erg/cm2 /s ' ^'^ ^ V 10^^ erg/s J \ 100 Mpc J ' ^^^^ 
where c^l is the luminosity distance to the source. This can be conipared with sensitivities of 
modern and planned 7-ray instruments. For example, iPaneque et al.l (|2007l ) gives the differential 
sensitivity of the GLAST LAT instrument, defined as 'the flux level over an energy interval of 
1/4 of a decade over which the statistical significance is 2 standard deviations in 1 year' roughly 
~ 0.8 X 10~^^ ergcm"^ s~^. This implies that GPS radio galaxies can be detected at 1 GeV photon 
energies by GLAST in its 1-year all-sky survey if dj^ < 100 (rje/VB)^^'^ L^^^^ LS^^J'^ ^gMpc. 
Thus, the most compact {LS < 100 pc) and powerful (Lj > lO^^ergs-^) sources with energy 
equipartition in their lobes fulfilled {rjc ^ ??b) are expected to be visible at GeV photon energies up 
to distances of ~ 1 Gpc. Note that even if the accretion disc luminosities in these objects are low, 
L\jy,4G < 1, the SSC component may provide a very strong source of target UV photons for the IC 
scattering. We also mention that since the sensitivities of EGRET and ^4 G/Lii' instruments are lower 
than the planned sensitivity of the GLASTLAT mission, one should not expect GPS radio galaxies 
to be detected by them. Obviously, there may be some exceptions, and it is possible that several 
unidentified high-latitude non- (or weakly) variable EGRET sources are associated with young 
radio objects. We emphasize that the estimates 27—28 provided above are very conservative, and 
any larger deviation from the energy equipartition, 77e 3> ?7b, as well as some other spectral shape of 
the electron injection function (e.g., the broken power-law discussed in the previous sections), will 
substantially increase the expected IC/UV flux. Interestingly, if the IC/UV component as discussed 
above (equation 28) extends to even higher-energy 7-rays, the expected IC/UV fluxes of GPS radio 
galaxies are accessible by the modern ground-based Cherenkov telescopes. For example, a rough 
detection limit for H.E.S.S. or MAGIC with a reasonable on-source exposure is ~ 10"^'^ ergcm^^ s^^ 
at 1 TeV photon energies, implying that GPS radio galaxies could, in principle, be detected by these 
instruments for source distances, again, dj^ < 100 {rjc/r]B)^^'^ L^^^^^ LS^qJ'^ L^y 4g Mpc. 



4. Conclusions 

There is an emerging agreement that compact radio sources characterized by inverted low- 
frequency radio spectra and symmetric ('lobe-lobe') radio morphology are young rather than frus- 
trated by an unusually dense ambient medium. That is because recent multiwavelength obser- 
vations reveal repeatedly that the properties of their environments and of their central engines 
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are very similar to the ones found in the extended radio sources. Here we utilize these newest 
multifrequency data, discussing evolution and non-thermal broad-band emission of the lobes in 
GHz-peaked-spectrum (GPS) radio galaxies. First, we propose a simple dynamical model for the 
GPS sources, exploiting the standard set of equations describing evolution of a relativistic jet — 
cocoon system. In contrast to the scenarios analyzed previously in the literature, however, we 
assume a uniform distribution for the ambient medium with number density no ~ 0.1 cm~^, as 
observed in central parts (< Ikpc) of elliptical hosts of radio galaxies. This gives a constant ad- 
vance velocity of the jet, which we fix as v^^ ^ 0.1 c, consistent with radio studies of the hotspots' 
proper motions in several GPS objects. With these few model assumptions we calculate all the 
other lobe parameters as functions solely of the jet kinetic power, Lj, and the source linear size, 
LS. In particular, the evaluated equipartition lobe magnetic field intensity i? ~ 1 mG is consistent 
with the value implied by the spectral-ageing analysis of radio emission from GPS galaxies, and the 
implied sideways expansion velocity of the lobes, Vc ^ 10^ cms^^, agrees with the outflow velocities 
of the line-emitting gas observed in many young radio objects. We note that both B and Vc depend 
weakly on the unknown jet power Lj, namely oc L^^^, which is a comfortable feature of the model. 

In the framework of the proposed dynamical description of young and compact radio sources, 
we follow the evolution of ultrarelativistic electrons injected from a terminal jet hotspot to the 
expanding lobes, taking the appropriate adiabatic and radiative energy losses into account. We 
find that, for a variety of the injected electron spectral shapes (which may be, in fact, differ- 
ent from the often assumed single power law oc 7~^), the resulting electron energy distribution 
within the lobes is expected to be of a broken power-law form, with critical electron energy 
7cr ~ 200 (Lj/10^^ ergs~^) . We argue that some plateau (or 'non-power-law curvature') around 
GHz frequency range is to be expected in the synchrotron spectra of GPS sources as a result of such 
a break. We also find that the total (unabsorbed) synchrotron luminosity of the lobes is expected 
to be constant in time during the GPS phase of the evolution (or, in other words, is expected to be 
independent on the source linear size), constituting about 10% (or more) of the total kinetic power 
of the jet. As for the formation of the inverted low- frequency radio spectra, we argue that it cannot 
be due to synchrotron self-absorption effects, because these cannot reproduce neither the slope nor 
the normalization of the observed fp oc LS~^'^^ anticorrelation. Instead, we favor free-free absorp- 
tion of radio photons by the neutral clouds of the interstellar medium, engulfed by the expanding 
lobes and photoionized by emission from the active center, as proposed previously by iBegelman 
(|l999l ). We note that these clouds may be naturally identified with the ones producing narrow- line 
emission (NLR) and also Hi absorption lines observed in many GPS radio sources. We speculate 
that both the observed anticorrelation of the neutral column density (corresponding to the observed 
Hi absorption lines) with the source linear size, A^hi oc LS~^''^^, as well as the anticorrelation of the 
spectral turnover frequency (being a result of a free-free absorption by the engulfed clouds) with 
the source linear size, is primarily due to the density decrease of the NLR clouds with the distance 
from the active center. Such a decrease, oc r~" with 1 < n < 2, is observed directly in several 
nearby Seyfert galaxies. 
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We estimate different photon fields within the lobes of GPS radio galaxies, due to emission from 
accretion disks, obscuring nuclear tori, synchrotron radiation of the lobe electrons, and starlight 
from the elliptical hosts. We calculate the resulting inverse-Compton components for different 
electron injection conditions, different jet powers, and different source linear sizes. We find complex 
high-energy spectra, extending from eV up to GeV (and possibly even TeV) photon energies. The 
resulting fluxes seem to be strong enough to dominate (or at least contribute significantly to) the 
radiation observed from GPS objects at UV and X-ray frequencies. We expect that in the case of 
GPS radio galaxies — where, due to large inclinations, the direct disk emission in the UV/X-ray 
range is expected to be obscured toward the line of sight by the dusty tori, while any jet radiative 
contribution is expected to be Doppler-hidden — the discussed non-thermal emission from the lobes 
is expected to be especially pronounced. In a subsequent paper, we investigate in more detail this 
issue through the analysis of the observed spectra of a sample of young and compact radio galaxies 
detected in the X-ray domain. Finally, we find that GPS radio galaxies are expected to be bright in 
the GeV photon energy range and that they can be detected by the forthcoming GLAST mission up 
to distances of the order of (conservatively) ~ 1 Gpc. Since the population of young radio sources 
is numerous when compared to the population of the extended ones, and since the high-energy lobe 
emission discussed in this paper is isotropic, one can expect that GPS radio galaxies are likely to 
constitute a numerous class of extragalactic GLAST sources. 

L.S. acknowledges support by the MEiN grant 1-P03D-003-29. L.O. gratefully acknowledges 
partial support from the INFN grant PD51, and thanks the Department of General Physics of 
the University of Torino for hospitality during the 2007 pre-fellow period. R.M. acknowledges 
support from MNiSW grant no. 1-P03D-009-28. J.K. acknowledges support by JSPS KAK- 
ENHI (19204017/14GS0211). L.O. and S.W. acknowledge support by BMBF/DLR through grant 
50OR0303. L.S., L.O. and S.W. were partly supported by the ENIGMA Network through the 
grant HPRN-CT-2002-00321. The authors thank G.V. Bicknell, M. Sikora, S.E. Healey, and the 
anonymous referee for very helpful comments. 



A. Electron Evolution 

Let us consider the kinetic equation 

= ^ {|7l A4(7, t)} + Q(7, t) (Al) 

describing the time evolution of the electron energy distribution A/'e(7, t) = A?e(7, t) V under the in- 
fluence of radiative (both synchrotron and inverse-Compton/Thomson-regime) as well as adiabatic 
energy losses, 

|7| = ^^7 + ciC/7', (A2) 
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where Q{'y,t) is the injection function, V is the volume of the system, V = dV/dt, c\ = AaT/SrUeC, 
and U = Ub + f^rad is the sum of the magnetic field and radiation field energy densities. Here we 
consider V = V{t), V = V{t), and U = U{t). 



Using the method of characteristics, one can find the solution to the kinetic equation Al in a 



form 



where 



and 7o = 'y{t = 0). Since 



lit) 



exp 








7o + ci J exp 


ft' 1 V 
J 3V 


dt" 


Udt' 



f * i ^ dt' 



Vo) 



where Vo = V{t = 0), one can find that A4 reads as 



7(t) = 



70 



-1/3 



l + ci7o/C/(^) 



-1/3 



Similarily, noting that 



,2ciJ*'YUdt' ^ 



1 + Cl 70 



dt' 



-1/3 



dt' 



one can rewrite A3 in the simpler form 



-A4(7,t) = 7-' (^) fQb{t')A lit'f {^^"^df, 



(A4) 



(A5) 



(A6) 



(A7) 



(A8) 



where we set He {lo-,t = 0) = 0. Note that in the above integral one has to substitute -^{t') with 
the expression A6, and then, after integrating, replace back 70 with 7. 
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Fig. 1. — Evolution of the electron energy distribution within the lobes of GPS sources for different 
jet powers (Lj = 10^^, 10^^, 10^^, and 10^^ erg s-^) and different source linear sizes (LS = 100 pc 
and 1 kpc; thick/upper and thin/lower lines, respectively). Two different source functions were 
considered for illustration, namely a single power-law Q{'j) oc 7"^-^ (dotted lines) or a broken 
power-law Q{'j) oc ^^^-^ for 7 < 7int and Q(7) oc 7^'^ for 7 > 7int (solid lines). In both cases. 



the minimum and maximum electron Lorcntz factors are 7n 



1 and 7n 



10^ while the 



normalization of the injection function is evaluated through the condition = r)eP with r]e = S 
and rjB = 0.3. Vertical dotted and dashed lines indicate critical electron energies 7cr = 200 rj^^ ^^45^ 
and 7int = mp/me, respectively. 
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Fig. 2. — Broad-band emission produced within the lobes of GPS sources with different jet kinetic 
power (Lj = lO^^ergs^^, lO^^ergs"^, lO^^ergs"^, lO'^^ergs"^) and different hnear sizes {LS = 
33 pc, 100 pc, and Ikpc). Illustrative parameters were considered: t/b = 0.3, ije = 3, Ly = 
1045 ergs-\ Luv = = 10^^ ergs'^ for Lj > 10^^ ergs-\ and Ljjy = Ljr = lO^^ ergs-^ for Lj < 
1045gpgg-i_ Single power-law injection function Q{'j) with spectral index s = 2.5 was assumed, 
with minimum and maximum electron Lorentz factors 7inin = 1 and 7max = 10^, respectively. 




Fig. 3. — Broad-band emission produced within the lobes of GPS sources with different jet kinetic 
power (Lj = 10^^ erg s^^, lO^^ergs"^, lO^^ergs"^, lO'^^ergs"^) and different hnear sizes {LS = 
33 pc, 100 pc, and Ikpc). Illustrative parameters were considered: t/b = 0.3, ije = 3, Ly = 
1045ergs-i, Luv = Lm = lO^^^ergs-^ for Lj > lO^^g^gg-i^ ^ ^^^^ ^ lO^^g^gg-i ^j. 

Lj < lO^^ergs"^. Broken-power-law injection function Q{'y) with si = 1.5, S2 = 3, break energy 
7int = rrip/me, 7min = 1, and 7i„ax = 10^ was assumed. 



